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Gluon excitations of the static quark potential and the hybrid 
quarkonium spectrum* 

K.J. Juge, J. Kuti, and C.J. Morningstar^ 

''Dept. of Physics, University of California at San Diego, La JoUa, California 92093-0319 

A comprehensive determination of the rich, low-lying spectrum of gluonic excitations in the presence of a static 
quark-antiquark pair is presented. Our results are obtained from several simulations on anisotropic lattices using 
an improved action and a large set of gluonic operators. The hybrid quarkonium states are studied in the quenched 
approximation using the Born-Oppenheimer expansion and nonrelativistic lattice QCD. 



1. Introduction 

Besides conventional hadrons, QCD suggests 
the existence of states containing gluonic excita- 
tions, such as glueballs and hybrid hadrons. Al- 
though conventional hadrons are reasonably well 
described by the constituent quark model, states 
with excited glue are still poorly understood. As 
experiments begin to focus on the search for glue- 
balls and hybrid mesons, a better understanding 
of these states from theory is needed. Due to the 
highly nonperturbative nature of the gluonic ex- 
citations in these states, lattice simulations offer 
at present the best means of theoretically probing 
glueballs and hybrid mesons. 

A great advantage in studying hybrid mesons 
comprised of heavy quarks is that such systems 
can be studied not only by direct numerical sim- 
ulation, but also using the Born-Oppenheimer 
(BO) expansion. In this approach, the hybrid me- 
son is treated analogous to a diatomic molecule: 
the slow heavy quarks correspond to the nuclei 
and the fast gluon field corresponds to the elec- 
trons||l|. First, one treats the quark Q and an- 
tiquark Q as spatially-fixed colour sources and 
determines the energy levels of the glue as a func- 
tion of the QQ separation r; each of these energy 
levels defines a potential ^^^(r) adiabatically. 
The quark motion is then restored by solving 
the Schrodinger equation in each of these static 
potentials. Conventional quarkonia arise from 
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the lowest-lying static potential; hybrid quarko- 
nium states emerge from the excited potentials. 
Once the static potentials have been determined 
(via lattice simulations), it is a simple matter 
to determine the complete conventional and hy- 
brid quarkonium spectrum in the leading Born- 
Oppenheimer (LBO) approximation. This is a 
distinct advantage over meson simulations which 
yield only the very lowest-lying states, often with 
large statistical uncertainties. 

Here, we present results for the spectrum of glu- 
onic excitations in the presence of a static quark- 
antiquark pair. Some of these potentials have 
been studied before This study is the first 
to comprehensively survey the spectrum in SU(3) 
gauge theory. Due to our use of anisotropic lat- 
tices and an improved action, we have been able 
to determine the static potentials for much larger 
values of QQ separation than previously studied. 
Using our potentials, we also determine the hy- 
brid quarkonium spectrum. Results from a pre- 
liminary nonrelativistic lattice QCD simulation 
are also presented. 

2. Computation of the potentials 

The first step in the Born-Oppenhcimer expan- 
sion is the determination of the rich spectrum 
of energy levels of the gluons in the presence 
of the quark and antiquark, fixed in space some 
distance r apart. At this point in the approxi- 
mation, the quark and antiquark simply act as 
static colour sources. The gluonic energies (or 
static potentials) may be labelled by the mag- 



Table 1 

Simulation parameters, including the coupling (3, 
input aspect ratio ^, approximate lattice spacing 
as, lattice size, and spatial link smearing param- 
eters (A and nx are defined in Ref. |^). 
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as (fm) 


Lattice 


(A,riA) 


2.2 


5 


0.27 


12^ X 48 


(0.10,4) 
(0.20,4) 
(0.30,4) 


2.4 


5 


0.23 


14^ X 56 


(0.10,8) 
(0.15,8) 
(0.20,8) 
(0.25,8) 


2.6 


3 


0.19 


10^ X 30 


(0.15,8) 
(0.30,8) 



nitude (denoted by A) of the projection of the 
total angular momentum of the gluons onto the 
molecular axis, by the sign of this projection (chi- 
rality or handedness), and by the behaviour un- 
der the combined operations of charge conjuga- 
tion and spatial inversion about the midpoint be- 
tween the quark and the antiquark. States with 
A = 0, 1, 2, . . . are typically denoted by the capi- 
tal Greek letters S, 11, A, . . ., respectively. States 
which are even (odd) under the above-mentioned 
parity-charge-conjugation operation are denoted 
by the subscripts g (u). The energy of the glu- 
ons is unaffected by reflections in a plane con- 
taining the molecular axis; since such a reflection 
interchanges states of opposite handedness, such 
states must necessarily be degenerate (A dou- 
bling). However, this doubling does not apply to 
the S states; E states which are even (odd) under 
a reflection in a plane containing the molecular 
axis are denoted by a superscript + (— ). Hence, 
the low- lying levels are labelled S+, S~, E+, 
Hg, Hu, Ag, Au, and so on. For convenience, we 
use r to denote these labels in general. 

Static potentials were extracted from Monte 
Carlo estimates of generalized Wilson loops. On 
a starting time slice the quark and antiquark 
were fixed at lattice sites a distance r apart. Sev- 
eral paths along the links of the lattice connect- 
ing the quark and the antiquark were then cho- 
sen, and our gluonic operators Of (to) were de- 




Figure 1. Examples of the paths from the quark 
to the antiquark used to construct the gauge field 
operators. 

fined as linear combinations of the path-ordered 
exponentials of the gauge field along these paths. 
Examples of these paths are shown in Fig. |[ The 
linear combinations were chosen such that the 
operators transformed irreducibly under all sym- 
metry operations. To reduce the mixings of our 
operators with excited states, the gluonic oper- 
ators were constructed from iteratively-smeared 
spatial links. We used the single-link smearing 
algorithm described in Ref. [^. The quark and 
antiquark then evolved in time, remaining fixed at 
their original spatial locations. To reduce statisti- 
cal noise, the static quark propagators, which are 
simply temporal Wilson lines, were constructed 
from thermally-averaged temporal links [Q, when- 
ever possible. The thermal averaging was done 
using the pseudo-heat-bath method (40 updates). 
At some final time slice to + r, evaluation of the 
gluonic operators O^^to+r) then completed the 
construction of the Wilson loops Wfj (r, r) . 

Monte Carlo estimates of the correlation matri- 
ces Wlj (r, t) were obtained in three simulations. 
In each simulation, several spatial link smearing 
schemes were used: one scheme was typically cho- 
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Figure 2. The static quark potential V^+{r) and 
some of its gluonic excitations in terms of the 
hadronic scale parameter tq against the quark- 
antiquark separation r. 



sen to work well for small r, another was opti- 
mized for large r. Various run parameters for the 
three simulations are given in Table |l|; other rel- 
evant parameters can be found in Ref. ||]. We 
used the improved action described in Ref. 
Our use of anisotropic lattices in which the tem- 
poral lattice spacing at was much smaller than 
the spatial spacing as was crucial. Configura- 
tion ensembles were generated using a mixture 
of Cabibbo-Marinari (CM) pseudo-heat-bath and 
SU(2) subgroup over-relaxation (OR) methods. 

The matrices W[j (r, r) were reduced in the data 
fitting phase to single correlators and 2x2 corre- 
lation matrices using the variational method de- 
scribed in Ref. § . These reduced correlators were 
fit using a single exponential and a sum of two 
exponentials in various ranges tmin to tamx of the 
source-sink separation. The two-exponential fits 
were used to check for consistency with the single- 
exponential fits, and in cases of favourable statis- 
tics, to extract the first-excited state energy. Best 
fit values were obtained using the correlated 
method. Error estimates were calculated using a 
1024-point bootstrap procedure. 
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Figure 3. The static quark potential V-^+ (r) and 
selected gluonic excitations (see Fig. ||). 

Our results for the static potential and its glu- 
onic excitations are shown in Figs. Results 
from the f3 = 2.2, P = 2.4, and (3 = 2.6 runs 
are shown using solid circles, squares, and trian- 
gles, respectively. The results are expressed in 
terms of the hadronic scale parameter ro- The 
definition of this parameter and a description of 
its calculation are given in Ref. The familiar 
static potential is shown as ; the solid curve is 
a fit to the data using a Coulomb plus linear form 

+ Sc/^ + • The curves for all other potentials 
are fits using cq -I- V^o + bir + b2r'^ ■ For all r stud- 
ied, the first-excited potential is the n„; hence, 
the lowest lying hybrid mesons must emerge from 
this potential. 

As r becomes very large, the linearly rising S+ 
potential suggests that the ground state of the 
glue may be modelled as a fluctuating tube or 
string of colour flux; in this picture, the gluonic 
excitations are expected to be phonon-like with 
energy gaps proportional to 1/r. However, it ap- 
pears that for r below about 1.5 fm, the gluonic 
spectrum cannot be explained in terms of a sim- 
ple string model. In Ref. jl], a QCD motivated 
bag model was successfully used to describe both 
the E+ and n„ potentials for a large range of r. 
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Figure 4. The static quark potential (r) and 
selected gluonic excitations (see Fig. |^). 

In this picture, the strong chromoelectric fields of 
the quark and antiquark repel the physical vac- 
uum (dual Meissner effect), creating a bubble in- 
side which perturbation theory is applicable. In 
the ground state, the inward pressure on the bub- 
ble from the physical vacuum balances the out- 
ward chromostatic force in such a way to produce 
a linearly confining potential. The addition of one 
or more gluons into the bag produces the excited 
potentials; the kinetic energy of the gluons inside 
the bubble is a key factor in determining the form 
of these potentials. This model has recently been 
revisited and results (in the ellipsoidal approxi- 
mation) for almost all of the potentials studied 
here are in remarkable agreement with our find- 
ings from the lattice simulations (see Ref. Q). 

3. Hybrid quarkonium 

The next step in the BO expansion is to re- 
store the quark motion by solving the radial 
Schrodinger equation. 



+ 2ii[E-V,s{r)] w(r) = 0, 



(1) 



where Vcff(r) = VQQ{r) + (L^-)/(2^r2), ^ is the 
reduced mass, and ^{r) — u{r)/r is the radial 



wavefunction. The total angular momentum of 
the meson is given by J — L -|- S, where S is 
the sum of the spins of the quark and antiquark, 
and the orbital factor L — ^qq + Jg , where 3g is 
the total angular momentum of the glue and ^qq 
is the orbital angular momentum of the quark 
and antiquark. In the LBO approximation, the 
eigenvalues L(L-l-l) and S{S+1) oiT? and are 
good quantum numbers. The centrifugal factor is 
then written as 



2A^ 



(2) 



For the S+ potential, (J^) = 0. For the n„ and 
S~ potentials, we assumed that the one gluon 
state was dominant with (J^) = 2. Mesonic 
eigenstates of parity and charge-conjugation are 
linear combinations of left- and right-handed glue 
states: |left) -I- e jright), where e = ±1. Let 
•q — ±1 denote the PC quantum number of the 
glue. Then in the LBO approximation, the parity 
(P) and charge conjugation (C) of each meson is 
given in terms of L and S according to 

P = £(-1)^+^+1, 

C = e 7? (-1)^+^+^. 



(3) 
(4) 



The potentials computed from Wilson loops in 
lattice simulations contain the self-energies of the 
temporal Wilson lines. These self-energy con- 
tributions are common to all of the static po- 
tentials and must be removed in order to ob- 
tain the VQ-Q{r) appearing in Eq. |^. Ideally, this 
can be done by measuring the S+ potential for 
very small r and comparing with the running 
Coulomb law as predicted from perturbation the- 
ory. In practice, this is difficult to do. Instead, 
we fit our results for the potential to the form 
VQ+ec/r-\-Kr; the constant Vq is then our estimate 
of the self-energy contributions to be removed. 

Results for the hh spectrum are shown in Fig. 
The scale was set using = 430 MeV as sug- 
gested from NRQCD simulations of hh and cc 
mesons (see Table XX of Ref. Q). The heavy 
quark mass Mjj was tuned in order to reproduce 
the experimentally-known T(15) mass: Mf = 
2Mb + Eq, where Eq is the energy of the lowest- 
lying state in the E+ potential. In the LBO ap- 
proximation, many mesons are degenerate: the 
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Figure 5. Spin-averaged bb spectrum in the lead- 
ing Born-Oppenheimer and quenched approxima- 
tions. Sohd hnes indicate experimental measure- 
ments. Short dashed lines indicate the S and 
P state masses obtained by solving the appro- 
priate Schrodinger equation in the Y,'^ potential 
using To 1 = 0.430 GeV and Mb = 4.60 GeV 
for the heavy quark mass. Long dashed and 
dashed-dotted lines indicate the hybrid quarko- 
nium states obtained from the n„ {L — 1,2) and 
E~ (L = 0, 1, 2) potentials, respectively. 



J-^*^ = '",1 S'-waves from the potential 
are degenerate; the 0"'""'",1"'""'",2"'""'",1"' P-waves 
from the E+ potential have equal masses; states 
such as 0-+,0+~,l++,l-+,l--,l+- from the n„ 
potential are also degenerate. 

Below the BB threshold, the LBO results are 
in very good agreement with the spin-averaged 
experimental measurements. Note that these re- 
sults make use of the quenched potentials (which 
ignore the light quarks) and do not include spin, 
retardation, and other relativistic effects. Above 
the threshold, agreement with experiment is lost, 
suggesting significant corrections from either the 
light quarks, relativistic effects, or possibly mix- 
ings between the states from the different adia- 
batic potentials. Note that the mass of the lowest- 
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Figure 6. Static potentials and radial probabil- 
ity densities against quark-antiquark separation 
r. The S+ static potential and the n„ and E" 
excitations are indicated by the dashed-dotted 
lines. The solid and short-dashed curves indi- 
cate the radial probability densities for the IS 
and IP states, respectively, corresponding to the 
results shown in Fig. |^. The extended nature of 
the lowest-lying hybrid state is shown by the 
radial probability density indicated by the long- 
dashed curve. 



lying hybrid (from the n„ potential) is about 
10.8 GeV. Hybrid mesons from all other static 
potentials are significantly higher lying. Above 
11 GeV, the LBO approximation based on the 
quenched Vqq potentials predicts a very dense 
population of hybrid states. The radial proba- 
bility densities for the 15 and IP conventional 
states are compared with that of the lowest-lying 
n„ hybrid state in Fig. ^ 

4. NRQCD simulations 

Hybrid quarkonium states may also be studied 
directly in numerical simulations. The nonrela- 
tivistic formulation of lattice QCD (NRQCD) §] 
is a particularly efficient means of carrying out 
such simulations. We have recently begun an in- 
vestigation of the hybrid quarkonium states using 
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Figure 7. Effective mass pfot sliowing tlie re- 
sults of a single-exponential fit to the correlation 
function of the hybrid T^~^ quarkonium state. 
The heavy quark propagates according to a spin- 
independent NRQCD action. 

a spin-independent version of the lattice NRQCD 
action. The action included only a covariant tem- 
poral derivative and the leading kinetic energy 
operator (with two other operators to remove 
0{at) and O(a^) errors); relativistic corrections 
depending on spin, E and B, and higher deriva- 
tives were not included. The action was chosen 
to correspond as closely as possible to the LBO 
approximation. In so doing, our treatments of the 
hybrid mesons using the Born-Oppenheimer ex- 
pansion and using NRQCD simulations differed 
in primarily two aspects: lattice spacing errors 
and retardation effects. The NRQCD simula- 
tions included retardation effects since the covari- 
ant Laplacian of the quark kinetic energy opera- 
tor was treated exactly. In contrast, the LBO 
approximation ignored all retardation effects by 
keeping only the leading term in the 1/c expan- 
sion of the covariant Laplacian. A comparison of 
results from the two approaches should afford a 
test of the adiabatic approximation. 

Because the signal-to-noise ratio was expected 
to be poor and the masses of the hybrid mesons 
were expected to be large, it was crucial to use an 



anisotropic lattice in which the temporal lattice 
spacing was much smaller than the spatial spac- 
ing. The operators used to calculate the static po- 
tentials were incorporated into our NRQCD sim- 
ulation code. We constructed large sets of gauge- 
invariant operators in order to minimize excited- 
state contamination of our hybrid meson correla- 
tors using variational techniques. We have com- 
pleted an initial run on a 12^ X 60 lattice using 
an aspect ratio ^ = 5 (that is, — 5at). As the 
purpose of this run was merely to discover the 
quality of signal which was possible, the quark 
mass was not tuned and a value a^Mb = 4.50 was 
used. The effective mass from this run for the 
1 ^ hybrid meson is shown in Fig. |^. A convinc- 
ing plateau was observed from 457 measurements, 
and we were able to extract the mass with a 3% 
statistical uncertainty. 

5. Conclusion 

A first comprehensive survey of the spectrum 
of quenched SU(3) gluonic excitations in the pres- 
ence of a static QQ pair was presented. The 
hybrid quarkonium states were calculated in the 
leading Born-Oppenheimer approximation. The 
effective mass for the 1 ^ hybrid from a prelim- 
inary NRQCD simulation was also shown. This 
work was supported by the U.S. DOE, Grant No. 
DE-FGG3-9GER4G546. 
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